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We study the statistics of curvature and torsion of Lagrangian trajectories from direct numer- 
ical simulations of homogeneous and isotropic turbulence (at 1Ze\ »s 280) in order to extract 
informations on the geometry of small scale coherent structures in turbulent flows. We find 
that, as previously observed by Braun et al |14| and by Xu et al [231 . the high curvature 
statistics is dominated by large scale flow reversals where the velocity magnitude assumes 
very low values. In order to focus on small-scales signatures, we introduce a cutoff on the 
velocity amplitude and we study the probability distribution of time— filtered curvature con- 
ditioned only on those events when the local velocity is not that small. In this way we are 
able to select small-scales turbulent features, connected to vortex filaments. We show that the 
conditional probability density of time— filtered curvature is well reproduced by a multifractal 
formalism. Finally, by studying the joint statistics of curvature and torsion we find further 
evidences that intense and persistent events are dominated by helical type trajectories. 
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1. Introduction 

The transport of particles in turbulent flows is an ubiquitous phenomenon rele- 
vant in very different physical contexts as cloud formation, pollutant dispersion, 
formation of "planetesimals" and turbulent mixing in chemical reactions [lr|4j] . To 
decribe it, the optimal framework results to be the Lagrangian one, where flow 
properties are studied following the trajectories of individual fluid elements. Thus, 
understanding the Lagrangian statistics of particles advected by a turbulent ve- 
locity field u(x, t) is a fundamental issue for both its theoretical implications (as, 
for instance, the development of stochastic models for dispersion and mixing) and 
for applications [l|, 0, 0]. Lagrangian description is also more suitable than the 
Eulerian one to reveal the geometry and the statistical "weight" of coherent struc- 
tures in turbulent flows. We cite here, for example, the preferential concentration 
of heavy/light particles inside hyperbolic/elliptic regions the strong ac- 

celeration bursts experienced by tracers when following streamlines around vortex 
filaments [13] or the regularising effects of inertia in the transition from viscous 
scales to inertial range scales in the behaviour of Lagrangian Structure Functions 
[l9l | . Studies of Lagrangian turbulence mainly involve the analysis of the statistics 
of velocity (or vorticity) and acceleration [1, 0, Ici-Hol. \v\ (for an updated recent 
review see Toschi & Bodenschatz [HI]). As suggested by Braun et al. [14], the tra- 
jectory of a Lagrangian particle x = x(i) can be also seen as a set of parametric 
equations representing a spatial curve in 3D. As such it can be described by means 
of tools borrowed from differential geometry: in particular, it is well known that 
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a curve is uniquely determined by the knowledge at every point (or equivalently, 
in our "time-like" parameterization, at every instant of time) of two fundamental 
geometric parameters: the curvature k and the torsion They are related to the 
velocity (u) and the acceleration (a) of the particle by: 

. , |u A a| aj_ . 

(where a± is the magnitude of the normal, or centripetal, acceleration) and 

„( t ) = !L££*>. (2) 
The torsion, unlike the curvature, has got a sign, which is related to the local 



helicity [12|, [13J. Since we are not dealing with helical flows (that is we have a 
symmetric distribution of regions with positive and negative helicity, hence the 
same is for torsion), from now on we will refer to torsion discussing its absolute 
value 9 = Dimensionally k and 8 are an inverse length, so we may expect that 
they assume very high values in correspondence of the small scale structures we 



are interested in, namely the vortex filaments. It has been shown [6|, [lfj, [17], |l_9(] 
that such filaments are of great relevance in the Lagrangian statistics of turbulent 
velocity and acceleration for both tracers and inertial particles (playing in the latter 
case a crucial role in the mechanism responsible for the strong inhomogeneous 
particle spatial distribution). The aim of this work is twofold. First, we will further 
characterize the geometry of these coherent structures by studying the statistics 
of extreme values of curvature and torsion from Direct Numerical Simulations, 
extending the previous analysis made in [14( to higher Reynolds numbers and 
supporting the numerical results with a phenomenological description in terms 
of the multifractal theory. Second, thanks to the high Lagrangian statistics, we 
will also present joint statistics of torsion and curvature, highlighting the intimate 
geometrical structure of small scales intense vorticity in turbulent flows. 

The paper is organized as follows. In section [2] we discuss the statistics of instan- 
taneous and time-averaged (and conditioned on velocity values) curvature, in the 
latter case deriving some analytical results in the framework of the multifractal 
formalism; in section [3] we analyze the joint statistics of curvature and torsion in 
order to stress its connection with the topology of small-scale coherent structures. 
Conclusions and perspectives are left to the last section. 



2. Curvature statistics 

We analyzed data produced in a Direct Numerical Simulation of homogeneous 
and isotropic turbulence in a cubic lattice of 1024 3 grid points [13], corresponding 
to 1Ze\ = 284, seeded with about two millions Lagrangian particles, evolving 
according to the dynamics 

^ = u(x(t),t), (3) 

where u(x, t), the fluid velocity field, is a solution of the incompressible Navier- 
Stokes equations 

d t u + (u • V)u = -VP + i/Au + f (4) 
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Figure 1. Log-log plot of the curvature PDF, normalized with the dissipative scale r). The solid line 
represents the k~ 5 / 2 scaling of the PDF for high curvature values. Data, here and hereafter, come from a 
high resolution DNS at 1024 3 collocation points, see text for details. 

V • u = 

which have been integrated on a triply periodic cubic box by means of a fully 
dealiased pseudospectral code (energy was injected at an average rate e by keeping 
the total energy constant in each of the first two wavenumber shells [2lf] ) . The 
equations ([3]) have been integrated linearly interpolating the values of the fluid 
velocity at the lattice sites. Let us stress that the acceleration field we can measure 
is continuous (and differentiable inside lattice mesh) because it is reconstructed 
from the local values of pressure gradients, laplacian of velocity and the external 
forcing along the trajectories. In such a way we can calculate a directly from the 
rhs of (|3|), without performing any differentiation of the particle velocity, thus 
minimizing noise effects, which are unavoidably present in experiments and, in 
general, in numerical simulations as well. 



2.1. Instantaneous and "time-filtered" curvature 

In figure (pQ) we show the PDF of the instantaneous curvature V(k). It is possible to 
notice the same power-law tail V(k) ~ k _5//2 , for large curvatures, observed in the 
numerical simulations by Braun et al. [14] and in the experiments by Xu et al. [23]. 
In the latter work it is argued that such a behaviour can be explained considering 
that the very high curvature events are linked to large scale flow reversals where 
the particle inverts its motion and its velocity goes to zero, u 2 ~ 0. Therefore, the 
PDF scaling for large curvatures is mainly determined by the Gaussian statistics 
of the velocity and, in this context, it can be reproduced analytically by simple 
arguments [23J. 

In order to avoid contamination from large scale motion, in [13] it was suggested to 
enhance the contribution of vortex filaments by studying the statistics of the cur- 
vature averaged over some time intervals A 22]]. The rationale at the basis of this 
operation consists in the observation that, while the "trapping" events of particles 
into the filaments are characterized by a certain degree of temporal persistency 
(due to the already stated spatio-temporal coherence of such structures), the flow 
reversals with very low velocity show very little autocorrelation times. Indeed, as 
also remarked in [231 ]. the larger is the curvature the lesser is the autocorrelation 
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Tabic 1. Conditional autocorrelation times and corresponding conditioning thresholds. The autocorrelation times 
have been extracted from the conditional autocorrelation function of the logarithm of curvature, defined as 

c , ) = (lQg(K(t + T))log(K(t))[log(«(t)) > logQO) - (log(K(t))|lQg( K (t)) > log(K*)) 2 -r/-r corr 
<log( K (t))2| log(K(t)) > log(**)) - <log( K (t))l log(«(t)) > log( K *)) 2 



At increasing the threshold (including in the computation only higher and higher curvature values) the autocor- 
relation decreases significantly. 
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Figure 2. Log-log plot of PDFs of instantaneous (•) and time averaged (x) curvature (over a window 
A = IOt,,). It must be noticed that only the tail at low k values has been effectively damped down by 
the filtering. Inset: Time series of instantaneous (solid line) and time— filtered (•) curvature relative to a 
trajectory in the high k tail of P(«(a=iQt ))• the arrow highlights the strong flow reversal event which is 
only very weakly damped by time averaging. 



time (in table Q] we report such times extracted from the autocorrelation func- 
tion of curvature, conditioned on its magnitude, together with the corresponding 
conditioning threshold). We have, thus, measured the PDFs of the time-filtered 
curvature: 

i rt+A/2 

«A(t) = -x / <t')dt' (5) 

A Jt-A/2 

for A = 7^,57^,107^, where t v is the Kolmogorov time scale. However, even for 
the greatest time "window" (10t„), we do not obtain the desired result. In fact, 
as it is evident from figure ([2]) where we report the PDFs of instantaneous and 
time-filtered (with A = 10r„) curvature, time-filtering is effective mostly at small 
k, where the PDF exhibits an important lowering of the tail. The explanation why 
filtering is not effective at high k values can be understood with the help of a 
single example. In the inset of figure (|2|) we analyse the effect of filtering around 
a very intense event in our statistics. We superimpose the temporal signal of the 
curvature together with the filtered one for a trajectory selected at random in the 
P(k(A=10t )) tail: at the instant of time at which the event of flow reversal occurs 
(indicated with an arrow) the curvature becomes so large that it is weakly damped 
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Figure 3. Log-log plot of conditional PDFs (of the time-averaged curvature), computed including only 
those events for which u 2 > A, for A = 0.02u^ ms (x), 0.1n^ ms (*) and 0.2u^ mB (□), respectively. The 
introduction of the cutoff induces a change in the slope of the far right tail (we also plot the unconditioned 
PDF (+) for comparison), reaching an asymptotic shape for the right tail insensitive to the increasing of 
A. 



even averaging over the longest time interval we used. 



2.2. Vortices identification and analytical results 

In order to really focus only on those high-curvature events given by large values 
of the numerator in (pQ) we propose to introduce also a cutoff (A) on the values 
of the velocity, i.e. we included in the computation of the filtered PDF only those 
events for which 



u 



> A. 



(6) 



We have tried various values A, see fig. ([3]), and observed that there is a tendency 
toward a PDF shape with a right tail insensitive to the increasing of the cutoff; 
keeping this in mind we will fix in what follows A = 0.02u^ ms . In this way, we 
define a conditional PDF: 



^a(^a) = f V(k a \u 2 > A). 



(7) 



As one can observe from figure (j4]), the application of this constraint results in a 
dramatic change in the shape of the PDF for high curvature values; in particular, 
the change of slope of Va(k(A=10t )) appears around the value 



KT] = K(A=10t„)»7 



0.1 



(8) 



where n is the Kolmogorov length. This new shape for the conditional PDF of the 
time-filtered curvature is compatible (once the two axes are rescaled accordingly) 
with the PDF of the time-filtered logarithm of the curvature (which was stud- 
ied in [13]): the reason of the convergence of the two procedures may be due to 
the fact that, by taking the logarithm of k, the two terms, the one coming from 
aj_ and the other coming from u 2 , contribute now linearly to the curvature, as 
log ft = log(aj_) — log(« 2 ); hence when one performs the time average (which is a 
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Figure 4. Log-log plot of the PDFs of instantaneous curvature without (+) and with (*) the cutoff on 
velocity values, and of the PDFs of time-filtered curvature without (x) and with (solid line) the cutoff. It 
is important to notice that the change of the PDF slope for large curvatures is due to the combination of 
the two operations (averaging in time and constraining the velocity magnitude). 



linear operation) the events at low u 2 values, characterized by short autocorrelation 
times, are selectively suppressed, similarly to what happens with the introduction 
of the cutoff. It seems then reasonable to infer that such change of the slope must 
be attributed to the presence of vortical structures, thus being connected to small 
scales turbulent properties. Indeed, because of the constraint imposed on u 2 , the 
highest curvature values are now due to large and temporally persistent accelera- 
tions (or, more correctly, to centripetal accelerations). If this is true, the large tail 
behaviour should be captured by the same multifractal arguments that have been 
proved successful for describing unconditional acceleration statistics pij . which is 
highly intermittent 0, d, 24]. We can write for the curvature PDF (denoted, from 
now on, with V(k), for the sake of simplicity) the following expression 



j j joint 



[a, u 2 )8 



4 ) <K« a 



A)dadu, 



(9) 



where 7 y j a int( a i u2 ) ls the joint PDF of a and u 2 , and the Heaviside function has 
been included to take into account the cutoff on the velocities. Note that we use 
a for aj_, neglecting geometrical constraints. Following the multifractal description 
presented in 24j we may estimate the joint PDF from the relation 



Vj i n t(a,u 2 ) = 8 (a - a(h,u))V(u 2 ); a(h,u) 



2h-l 3 
V 1+h U X + h Lr 



3h 
1+h 



(10) 



where v is the viscosity, Lq is the integral scale (we will keep Lq ~ 1 in the 
following) and h the scaling exponent characterising the local velocity fluctuation. 
Let us stress that in the expression of a(h, u) a possible over-all stochastic prefactor, 
Reynolds independent, can be included, in order to change the degree of correlation 
between the acceleration and the local velocity field. Such a prefactor would change 
a bit the joint PDF, but not the global shape of the result. At the Kolmogorov 
scale, the probability to have an exponent h is given by [U 



(11) 
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where Tl is the Lagrangian integral time scale and D{h) is the fractal dimension 
of the set characterised by the exponent h. Using for the probability density of u 2 
a third order \ 2 distribution: 

V(u 2 ) oc u 2 e~ u2 / 2 (12) 

we obtain: 



/a n 3 — D(fl) ^ 

da du dh(pj 1+h u 2 e~^ 5{a - u^)5 [k - -^j @{u 2 - k). 

(13) 

Integration over a and u yields the final expression 

h / 9 2(l + h) \ 

. . I nma * ( 1 + h \ n 5h-l + D(h) / //J 1-21. \ 

(k) A \T^zh) l ~ 2h exp [ 2 — ) ( } 

valid for n ~ 2.0<r K (correspondingly kt) ~ 0.03)0 Working out the integral with 
h € [0.16,0.38] (for integrability reasons exposed in (24[) and using for D{h) the 
empirical form proposed by She and Leveque [2^] (which was also shown to fit very 
well experimental data [26J], at least in the range of singularity exponents h that 
we are spanning): 

p(h) 

D{h) = l+p{h)(h-^j+2{^j 3 (15) 

(where p(h) = (3/ ln(2/3)) ln[(9/i — l)/61n(3/2)]), we have the result shown in 
figure ©, in excellent agreement with the measured PDF for large curvatures. 
More in detail, it is worth remarking that the two curves are almost coincident for 
values of curvature larger than the "critical" value R. Let us notice that the good 
agreement between data and the multifractal formalism only after the introduction 
of a cutoff in the single point velocity is not an accident. We must expect that 
phenomenological models based on the energy cascade mechanism work well only 
for those events that are not dominated by smooth quasi-laminar fluctuations, as 
it would be the case for the whole curvature statistics without cutoff. 



Joint statistics of curvature and torsion 



Torsion measurements are very difficult in experiments 23 1 and available only 
at very low Reynolds numbers in previous numerical works 14|. The difficulty is 
due to the fact that they involve the time derivative of the acceleration. In our 
numerical data-base acceleration is differentiable inside mesh grids and continuous 
everywhere: we could therefore measure torsion with a good signal to noise ratio. 

It is natural to suppose that the emergence of coherent structures with peculiar 
geometries induces somehow a correlation between curvature and torsion of tra- 
jectories, as they contain informations on such structures. We studied, then, the 



joint statistics of the two quantities. We calculated the residual value [22| of the 



lr This domain limitation arises from the calculation done with the function. 



September 11, 2009 16:12 Journal of Turbulence scagliarini'curvature'arXiv 



10 



Taylor & Francis and I. T. Consultant 



a 

Q. 



10" 



10" 



10"' 



10"' 



10" 



10"' 



DNS • 

multifractal prediction 



0.05 




0.45 



Figure 5. Log-lin plot of the PDF with cutoff on u' 2 of the time filtered curvature (A = IOt,,) and the 
corresponding multifractal prediction. 
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log(Kti) 

(a) 

Figure 6. Residuals of the logarithms of curvature and torsion, instantaneous (a) and time-averaged (b) 
over a time window A = IOt,,, plotted as a 2D map with isocontours (only TZ > is shown). Correlation 
is lesser on exterior contours. In the right panel we also plot a straight line parallel to the plane bisector 
log(0) = log( K ). 



logarithms of k and 9, that is the quantity 

ft«log(«)) A , <log(0)) A ) = ?W(log(«))A, aog(0))A)-P(<log(K))A)Wog(0)>A) 

(16) 

which quantifies the degree of correlation existing between the two variables. 
In figure © we plot ft(log(/c), log(0)) and ft((log(K)) A , (log(<9)) A ), with A = 
10t v , projected onto the log(0) vs log(«;) plane (isocontours are drawn for pos- 
itive correlation 1Z > 0). In figure (|7|) we also show the original joint PDF 
*Pjoint({^og(K))/\, (log(6 , )) A ) for both the instantaneous and time-filtered quanti- 
ties. The most relevant feature to be stressed is that the effect of time filtering is 
to accumulate the isocontours to intensify the correlation along the straight lines 
log(0) = log(«;) + c (parallel to the plane bisector). This equation is exactly what 
one obtains by taking the logarithm of 



(17) 
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Figure 7. Joint probability distribution functions of the logarithms of curvature and torsion, instantaneous 
(a) and time— averaged (b) over a time window A = IOt,,. 

which is the relation linking torsion and curvature for a "perfect" cylindrical helix 
(A being its reduced step and R the cylinder radius), that is the prototype of 
spatial curve that we imagine when thinking of a vortex tube. From the value of 
the intercept we can then also estimate the ratio between the step and the radius 
of the "mean" helix to be A/r ~ 3. These two geometrical parameters are related 
to the inverse curvature (i.e. the radius of curvature) via 



r 2 + A 2 



(18) 



furthermore from the critical value at which the change of slope of the curvature 
PDF occurs (discussed in the previous section), we get for the radius of curvature 
R ~ IO77. Hence, substituting the measured ratio A/r into (I18j) . we obtain, as an 
estimate for the radius of the helix, 



r ~ r\. 



(19) 



These observations strengthen furthermore the conjecture that when filtering and 
avoiding the very low velocity events, the study of the statistics of curvature and 
torsion provides an insight on the signatures of vortex filaments. 



4. Conclusions 

We have analyzed data from DNS of Lagrangian tracers, focusing our attention on 
the statistics of parameters able to characterize the geometry of particle trajecto- 
ries, namely curvature and torsion. We found that the statistics of instantaneous 
curvature is mainly dominated by large scale flow reversal events, rather tha n by 
small scale structures like vortex tubes as also found in experimental works in [23j. 
Indeed, the asymptotics of the PDF can be reproduced by simple Gaussian statis- 
tics arguments [231 ] . In order to unravel small-scales turbulence features, we meant 
to compute the PDF of curvature averaged in time, including only those events for 
which the squared velocity magnitude u 2 was greater than a certain threshold, so to 
avoid the reversals (where u 2 ~ 0) and to exalt the statistical signature of vortices. 
This approach successfully resulted in a marked change of slope in the PDF tail at 
high curvature values, pointing out the possible presence of vortex filaments. The 
criterion here presented is therefore different from what proposed in other works 
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15J ]. Moreover, we could derive an analytical expression for Va(ka)j m the context 



of the multifractal formalism, which fits well the far tails, providing an interesting 
matching between geometry and phenomenology. The analysis of joint statistics of 
curvature and torsion has been also addressed, confirming our statements on the 
topology of coherent small scale structures. 
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